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Abstract 

Using a sample of 227 million T(4S) — > BB events collected with the BABAR detector at the PEP-II 
B Factory in 1999-2004, we study B~ -> D° K* (892) " decays where K*' -> K Q S n~ and D° -> 
K~ vr+, K~ vr+ vr°, K~ 7T+ 7T+ tt~ (non-CP final states), K+K~, vr+ vr" (CP+ eigenstates), K° vr°, 
Kg and to (CP— eigenstates). The partial rate charge asymmetries Acp and the ratios TZcp 
defined in the literature as the sum of the B + and B~ partial rates to a charged K* and a D° 
CP-eigenstate divided by the B — > D° K* decay rate, are sensitive to the angle 7 of the CKM 
unitarity triangle. We measure: 

A C p+ = -0.09±0.20(stat.)±0.06(syst.) 

A C p- = -0.33 ± 0.34(stat.) ± 0.10(syst.) (+0.15 ± 0.10) • (A C p- - A C p+) 

n C p+ = +1.77 ± 0.37(stat.) ± 0.12(syst.) 

K C p- = +0.76 ± 0.29(stat.) ± 0.06(syst.) 



0.04 
0.14 



The third uncertainty quoted for the CP— measurements reflects possible interference effects in the 
final states with cf) and uj resonances. All results are preliminary. 
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1 INTRODUCTION 



The measurement of CP violation in B meson decays offers the means to over-constrain the unitarity 
triangle. A theoretically clean measurement of the angle 7=arg(— V u dV* b /V c dV* b ) is provided by 
the B — ► Z)(*)° [1] decay channels in which the favoured b — ► cms and suppressed b — > ucs 

penguin- less processes interfere [2]. In this paper we implement the method proposed by Gronau, 
Wyler and London by looking at the interference between B~ — ► D° K*~ and B~ — > D° X* - when 
both B>° and B>° decay to a CP eigenstate. 
We define [3]: 

r(g- -» z%p ± g*-) - r(B+ -» z^ P± if *+) 
CP± r(B- -► d° cp± k*-) + r(5+ £>£ p± jf*+) ' > 
T(B- D° CP± K*-) + T(B+ -> D cp ± K*+) 

n cp± = t(b~TWk^) • (2) 

Acp± are the direct CP-violating asymmetries associated to each D° CP eigenstate. 1Zcp± is twice 
the ratio of the charge-averaged branching fraction of the B — > D°cp± K* decays to that of B~ 
— > D° if* - . The latter branching fraction has been measured previously by BAR4R [4], CLEO [5] 
and Belle [6]. These measurements average to (5.7 ± 0.6) • 10 -4 . Both Acp± and TZqp± carry 
CP-violating information. Neglecting D°-D° mixing: 

1Zcp± = 1 =t 2rs cos 5 cos 7 + r 2 B (3) 



Acp± x TIcp± = ± 2rg sin <5 sin 7 (4) 

where (5 is the CP-conserving strong phase phase difference between the B — > D° X* - (favoured 
B decay) and B~ — > B>° if*" (suppressed) amplitudes, 7 is the CP-violating weak phase difference 
in the same processes and re ~ 0.1-0.3 is the magnitude of the ratio of the amplitudes [3]. 

2 THE BABAR DETECTOR AND DATASET 

We look for B~ — > D^p If *~ decays with the data collected near the T(4S) resonance with the 
BABAR detector at the PEP-II storage ring between October 1999 and May 2004. We accumulated 
an integrated luminosity of 205 fb _1 on the peak (227 million BB pairs) and 16 fb _1 40 MeV below 
the resonance. The BABAR detector is described elsewhere [7]. We focus on the components which 
are relevant to this analysis. Charged-particle trajectories are measured by a five-layer double sided 
silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) immersed in a 1.5 T solenoidal 
magnetic field. Charged-particle identification is achieved by combining the light measured in a 
ring- imaging Cherenkov device (DIRC) with the ionization energy loss (dE/dx) measured in the 
DCH and SVT. Photons are detected in a CsI(Tl) electromagnetic calorimeter (EMC) inside the 
coil. We use GEANT4 [8] based software to simulate the detector response and account for the 
varying beam and environmental conditions. 



3 ANALYSIS METHOD 



We reconstruct B~ -► D° K* (892)" decays. We select K*~ candidates in the K*~ -► if° tt~ 
mode and D° candidates in 8 decay channels. We optimize our event-selection criteria to minimize 
the statistical error on the signal yield, determined for each channel using simulated signal and 
background events. In order to maximize statistics and as we measure ratios (asymmetries and 
ratios of branching fractions), we use rather relaxed criteria to define a track or a photon candidate. 
In general a track must originate from the interaction point within 1.5 cm in the transverse plane 
and 10 cm along the beam. Very soft tracks, measured only by the SVT are also used. Loose 
criteria are employed to separate charged pions and kaons. Pions (kaons) are identified with more 
than 95% (90%) efficiency by an algorithm which rejects more than 85% (95%) of the kaons (pions) 
in the geometrical acceptance of the DIRC. Photon candidates are isolated energy deposits in the 
EMC with an energy above 30 MeV and a shape consistent with that of a photon-induced shower. 

Kg candidates are made from oppositely charged tracks assumed to be pions. K® used to 
make D°cp± (K*) candidates are mass (mass and vertex) constrained. The pion pairs are selected 
according to their reconstructed invariant mass, required to be within 13 MeV/c 2 from the known 
value [9]. For those used to search for a K* we further require their flight direction and length 
to be as expected for a Kg coming from the interaction point. The Kg candidate flight path and 
momentum have to make an acute angle. The Kg candidate flight length in the transverse plane 
has to exceed its uncertainty by 3 standard deviations. 7r° candidates combine pairs of photons 
with a total energy greater than 200 MeV and an invariant mass between 115 and 150 MeV/c 2 . A 
mass constraint fit is applied to the selected 7r° candidates. K* candidates combine a Kg and a 
charged particle which are constrained to come from a single decay vertex. We select K* candidates 
which have an invariant mass within 75 MeV/c 2 of the known value [9]. Finally, since the K* in 
B —>■ D° K* is polarized, we require |cos 9x*hel I > 0.35, where 0K*hei is the angle in the K* rest 
frame between the daughter pion and the parent B momenta. Here and in other occurrences the 
helicity distribution discriminates well between a B meson decay and an event from the e + e~ — ► 
qq continuum. 

We select D° in even CP (CP+) eigenstates (K + K~ , ir + n~), odd CP (CP-) eigenstates 
(Kg 7T°, Kg (f) and Kg lo), and non-CP (flavour) eigenstates (K~ ir + , K~ tt + tt° and K~ ir + ir + 
7r - ). Composite particles included in the CP— modes are vertex constrained, (p (lo) candidates 
are constructed from K + K~ (ir + ir~ ir°) particle combinations with an invariant mass required 
to be within 12 (20) MeV/c 2 of the known values [9]. We further select the lo candidates with 
requirements on two helicity angles. The normal helicity angle, between the D° momentum in 
the rest frame of the lo and the normal to the plane containing all 3 decay pions, must have its 
cosine above 0.25 in magnitude. The Dalitz-helicity angle [10], defined as the angle between the 
momentum of one daughter pion in the lo rest frame and the direction of one of the other two pions 
in the rest frame of those two pions, must have a cosine with a magnitude less than 0.9. 

Except for the Kg 7r° final state, all D° candidates are mass and vertex constrained. We select 
D° candidates with an invariant mass differing from the known mass by less than 12 MeV/c 2 for 
all channels except K° s vr° (30 MeV/c 2 ) and K° lo (20 MeV/c 2 ). These limits are about 2 standard 
deviations of the nearly-Gaussian resolutions. 

To suppress the background due to e + e~ — > qq reactions, q G {u,d,s,c}, we require |cos 6b\ < 
0.9 as, in the T(4S) rest frame, B mesons are produced with a sin 2 0b distribution of the angle 9b 



of their momentum with the beam axis. We also use global event shape variables which translate 
quantitatively the fact that in the T(4S) rest frame, qq continuum events have a two-jet like 
topology whereas BB events are more spherical. We require |cos 6thrust\ < 0.9 where Othrust is 
the angle between the thrust axes of the B candidate and the rest of the event (i.e. what is left 
after removing the tracks and clusters associated with the B candidate). We construct a Fisher 
discriminant [11] from cos 6 thrust and Legendre monomials [12] describing the energy flow in the 
rest of the event. 

We identify final B candidates using two [nearly] independant kinematic variables: the beam- 
energy-substituted mass tties = \J (s/2 + po • pb) 2 /-^o ~~ P% an< ^ tne ener gy difference AE = 
E* B — y/s/2, where the subscripts and B refer to the e + e~-beam-system and the B candidate 
respectively; s is the square of the center-of-mass (CM) energy and the asterisk labels the CM 
frame. For all signal modes, the ttt-es distributions are described by the same Gaussian function 
Q centered at the B mass with a 2.6 MeV/c 2 resolution (except tt° which has a 2.8 MeV/c 2 
resolution). The AE distribution depends more strongly on the signal mode. AE is centered on 
zero for signal with a resolution of 11-13 MeV for all channels except Kg ir° for which the resolution 
is asymmetric and about 30 MeV. We define a signal region \AE\ < 25 MeV for all modes except 
K® 7T° which uses < 50 MeV. 

For those events consistent with more than one B candidate — this occurs in < 25% of selected 
events depending on the D° mode — we choose that with the smallest \ 2 formed from the differences 
of the measured and true D° and K* masses scaled by the mass spread which includes the resolution 
and, for the K*, the natural width. Simulations prove that no bias is introduced by this choice and 
the correct candidate is picked at least 82% of the time. 

The total reconstruction efficiencies 6 after corrections, according to simulation of signal events, 
are: 12.5% and 12.9% for the CP+ modes D° -► K + K~ and vr + vr"; 5.0%, 10.2% and 2.3% for the 
CP- modes D° -► K® vr°, K° s <f> and K° s uj; 13.5%, 5.2% and 8.0% for the D° -► K~ vr+, K~ vr+ 
7T° and K~ ir + ir + ir~ non-CP modes. 

4 PEAKING BACKGROUNDS 

To study BB backgrounds we look in sideband regions away from the signal region. We define a 
AE sideband in the interval -100 < AE < -60 and 60 < AE < 200 MeV for all modes except D° 
— > Kg 7T° for which the inside limit is ±95 rather than 60 MeV. The lower limit (—100 MeV/c 2 ) 
is chosen to avoid selecting much of the background coming from D*° — > D° 7r°/7 decays. In this 
AE sideband we see no significant evidence of a peaking background which could leak into the 
signal region. We define a second control region in the m D o sideband away from the D° mass peak. 
This provides sensitivity to background sources which mimic signal both in AE and tties: the 
(doubly) peaking background. This pollution comes from either charmed or charmless B meson 
decays which do not contain a true D°. As many of the possible contributions to this background 
are not well known, we attempt to measure this contribution by including the m D o sideband in the 
fit. Relevant plots are shown on the third row of Figure 1. 

A notable background for the B~ — ► (D° — > 7r + 7r~) {K*~ — > K® ir~) mode is the dacay 
B~ — ► (Kg7r + Tr~) D oTr~~ which contains the same final state as the signal but has a 600 times higher 

6 The branching fractions of the unstable daughter particles are not included. 



branching fraction. We therefore explicitly veto any selected P-candidate containing a (K® ir + tt ) 
combination within 25 MeV/c 2 of the D° mass. 

5 YIELD AND CP FIT 

An unbinned extended maximum likelihood fit to tttes distributions in the range 5.2 < tttes < 
5.3 GeV/c 2 is used to determine yields and CP-violating quantities Acp and IZcp- We use a uni- 
versal Gaussian function Q to describe the signal shape for all modes considered. The combinatorial 
background in the ttt-es distribution is modeled using a threshold function Arg first considered by 
the Argus collaboration [13]. Its shape is governed by one parameter £ that is left to float in the 
fit and a second, Emax = y/s/2, which is fixed at 5.2910 GeV/c 2 . We fit simultaneously tttes 
distributions of nine samples (1) the non-CP, (2) CP+ and (3) CP— D° modes in the (i) signal 
region, (ii) the m D o sideband and (Hi) the AE sideband. We fit three types of probability density 
functions (PDF) combining (tttes; Xj) and Arg(rriEs', weighted by the unknown event yields. 
We call Xj, the mean and standard deviation of Q and £ the shape parameters of the PDF . In the 
AE sideband, we fit: Arg. In the m^o sideband we fit: a s b • Arg + b s b • Q. where Q accounts for 
the doubly peaking f?-decays. Finally, we fit: a ■ Arg + b ■ Q + c • Q in the signal region. Here b 
is scaled from b s b according to the m D o sideband and signal window widths and is fixed; c is the 

— > D° K*^ signal. The motivation to perform the fit just described is as follows: the non-CP 
mode sample with relatively high statistics helps constrain the PDF shape for the low statistics 
CP mode distributions. The AE sideband sample helps define the Argus background shape. 

The values of £ obtained for each data sample were found to be consistent with each other, 
but subject to large uncertainties. We have therefore constrained £ to take the same value for all 
data samples in the fit. We assume that the B decays found in the m D o sideband have the same 
final states as the signal and that we can therefore fit the same Gaussian to the signal and doubly 
peaking B background. 

Furthermore, the CP samples (in the signal region) are split into the B~ and B + sub-samples. 
The likelihood function is written so as to directly extract the four physical quantities, A+, A-, 1Z+ 
and 7Y_. The doubly peaking P-background is assumed to not violate CP and so is split equally 
between the B~ and B + sub-samples. This assumption is considered further when we discuss the 
systematic uncertainties. Fig. 1 shows graphically the results of the fit. The nominal fit results are 
shown in table 1. 





non-CP 


CP+ 


CP- 


Yield 


498 ± 29 


34.4 ± 6.9 


15.1 ± 5.8 


Npeak 


10.9 ± 6.6 


2.5 ± 1.3 


2.4 ± 2.2 


A C p 




-0.09± 0.20 


-0.33± 0.34 


Hop 




1.77± 0.37 


0.76 ± 0.29 



Table 1: Results from the fits. For each D mode class, we give the event yield, the peaking 
background contribution, Acp and IZcp- The uncertainties are statistical only. 

The fit stability and accuracy have been studied by performing 1000 simulated experiments at 
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Figure 1: The simultaneous fit. The tties distributions are shown for the non-CP modes in the first 
column, CP+ modes in the second column and CP— in the third. The fits in the signal region are 
shown in the top two rows, fits of the peaking background in the m^o sideband are in the third 
row and the fits of the background in the AE sideband are in the fourth. For the CP modes in the 
signal region, B + and B~ distributions are shown separately in the first and second rows. 



the optimum found by the fit in the parameter space. The value of the likelihood function in the 
fit falls well within the range observed with the simulated experiments. 

The yields we observe for the non-CP modes are compatible with our previous measurements 
of the B~ — > D° K*~ branching ratio [4]. The statistical significance of the CP+ (CP—) yields are 
6.8 (2.9) standard deviations. It should be noted that the measured Acp- could differ from the 
physical quantity, as background under the 4> and lu resonances with odd and even CP could be 
present. These possible interference effects are accounted for in the systematics. 

6 SYSTEMATIC STUDIES 

6.1 CP-asymmetries 

Although most systematic errors cancel in the asymmetries, three effects must be considered. An 
asymmetry inherent to the detector or analysis may exist. After running the analysis code on a high 
statistics B~ — > D° ir~ sample (with K* cuts removed), the final selection shows an asymmetry of 
(—1.9 ± 0.8)%. We quote a systematic uncertainty of ± 2.7%. 

The second substantial systematic effect is that due to a possible CP-asymmetry in the peaking 
background. Although there is no physics ansatz that requires the peaking background to be 
asymmetric, it cannot be excluded. To get an estimate, we remark that if there is an asymmetry 
■Apeak-, a systematic error on Acp is given by A pea k x \-, where b is the contribution of the peaking 
background and c the signal yield. Assuming conservatively A pea k < 50 %, we obtain systematic 
errors of ± 3.5 % and ± 8.0 % on Acp+ and Acp- respectively. 

Finally a systematic correction must be applied to Acp- which stems from the background under 
the 4> and oj resonances. Looking at the D° — > K~ K + Dalitz plot [14] we see a contribution of 
the ao(980) under the <fi. The ao(980) is a ++ state which induces a CP+ pollution. We crudely 
estimate a bias on Acp- with a linear dependance on the difference of the measured asymmetries, 
8A_a cp _ = 0.15 ± 0.10 • (Acp- — Acp+) which we quote as a separate systematic uncertainty. 

6.2 n CP 

IZcp are ratios of rates of processes differing by the final state of the D°. We have to consider 
the uncertainties affecting the selection algorithms for the differing D channels. They bring small 
correction factors which account for the difference between the actual detector response and the 
simulation model. The main effects stem from the approximate modeling of the tracking efficiency 
(1.2 % per track), the K® reconstruction efficiency for CP— modes of the D° (2.0 % per K®), the 
7T° reconstruction efficiency for the K® 7r° channel (3 %) and the efficiency and misidentification 
probabilities from the particle identification (2 % per track). A substantial effect is the uncertainty 
on the measured branching ratios [9]. We find 4.8 % and 7.5 % for the systematic uncertainties 
on the selection of the D° to CP+ and CP— channels and 4.0 % for D° to non-CP channels. 
Altogether, we find systematic uncertainties on the relative efficiencies to be 6.2 % and 8.5 % on 
TZcp+ and IZcp- respectively. Here also we quote as a separate systematic uncertainty, a possible 
bias due to structure under the 4> and to resonances. We conservatively estimate this bias to lie 
between -18 and -5 % using the D° -► K~ K + K Q S Dalitz plot [14]. 



7 PHYSICS RESULTS 



We quote the final results: 

A CP+ = -0.09 ± 0.20(stat.) ± 0.06(syst.) 

A C p- = -0.33 ± 0M(stat.) ± 0.10(syst.) (+0.15 ± 0.10) • (A C p- - A C p+) (bias) 

n C p+ = +1.77±0.37(s*at.) ±0.12(syst.) 

^ CP _ = +0.76 ± 0.29(s*at.) ± 0.06(syst.) I (6ms) 

For the CP— measurements, the third uncertainty reflects a possible bias due to interference 
effects in the final states with <fi and u resonances. 

These can be compared with the preliminary results by the Belle collaboration [6] using 95.8 
BB pairs: A C p+ = -0.02 ± 0.33 ± 0.07; A C p- = +0.19 ± 0.50 ± 0.04. The results are in agreement 
at the one standard deviation level. As expected our statistical errors are reduced. We quote 
much larger systematic uncertainties because we account for a possibly high CP-asymmetry in the 
background and for structures under the (j) and to resonances whereas reference [6] does not. 

We use equation (3) to derive r 2 B = 0.23 + 0.24. Although the central value is large (it translates 
to rs=0.47), a null value is within one-standard deviation. 

8 SUMMARY 

In summary, we present preliminary observations of the decays of charged B mesons to a K* and 
a D° where the latter particle is seen in final states of even and odd CP. We express the results 
in terms of V-cp and Acp. With more statistics, these quantities will constrain r#, the ratio of 
amplitudes defined in equations (3 and 4) and the angle 7 of the unitary triangle by application of 
the Gronau-London-Wyler method. 
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